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1.  Introduction 


In  this  work,  radial  wavefunctions  describing  the  quantum  behavior  of  a  spin  one-half  particle  far 
from  the  center  of  a  uniform  atomic  waveguide  are  formulated.  The  model  waveguide  potential 
is  an  ideal  transverse  quadrupole  magnetic  field  that  extends  to  infinite  distances  and  is 
unbounded  in  magnitude.  The  spin  of  the  particle  can  be  considered  in  either  a  global  basis  or  in 
a  local  basis.  The  quantization  axis  for  the  global  basis  is  directed  in  the  z-direction  along  the 
axis  of  the  guide  while  the  quantization  axis  for  the  local  basis  is  taken  to  be  directed  everywhere 
along  the  local  direction  of  the  quadrupole  magnetic  field.  In  either  basis,  the  atomic 
wavefunction  has  two  components  one  of  which  might  seem  more  “trapped”  than  the  other. 
However,  in  each  of  these  bases  there  is  a  significant  component  of  the  spinor  eigenstates  that 
extends  to  infinity  and  represents  the  amplitude  for  an  atom  to  be  found  outside  of  the  trapping 
region  of  the  waveguide. 

The  asymptotic  solutions  developed  here  describe  the  behavior  of  the  radial  wavefunctions  at 
distances  far  from  the  guide  center  and  thus  provide  a  means  of  extending  the  series  solutions 
derived  in  Golding  (2010)  to  these  regions.  The  combination  of  the  two  series  can  be  combined 
into  an  efficient  and  complete  solution  to  the  waveguide  eigenstate  problem. 

In  the  global  spin  basis,  both  spinor  components  of  the  radial  solution  appear  to  be  unbound. 

They  oscillate  and  only  gradually  decay  along  the  radial  direction  outward  towards  infinity. 
However,  in  the  local  spin  basis  one  component  of  the  spinor  appears  to  be  bound  and  the  other 
component  appears  unbound.  This  results  from  a  nearly  exact  cancellation  of  the  global 
components  far  from  the  center  when  viewed  in  the  local  basis.  This  behavior  is  clearly  related  to 
the  adiabatic  approaches  often  used  to  describe  these  waveguides  when  the  axial  bias  field  is 
made  large  (Sukumar,  1997;  Brink  2007).  In  the  adiabatic  approach,  the  unbound  component  is 
effectively  considered  to  be  uncoupled  from  the  bound  component  and  ignored  in  a  first 
approximation. 

A  goal  of  our  research  is  to  precisely  understand  the  behavior  of  atoms  guided  on  atom  chips  and 
moving  through  atom  chip  devices.  A  complete  description  of  the  solution  to  the  waveguide 
problem  is  a  required  starting  point  in  the  pursuit  of  this  goal. 

This  report  begins  by  laying  out  the  radial  equations  that  must  be  solved  to  describe  the 
transverse  modes  of  an  atom  guide.  The  leading  behavior  of  the  solutions  to  these  equations  at 
large  distances  is  then  developed.  To  properly  connect  the  four  solutions  that  represent  the 
leading  behavior  the  phase  integral  technique  is  used.  The  phase  integral  technique  essentially 
produces  the  relative  phase  and  amplitudes  needed  to  create  an  asymptotic  expansion  that 
correctly  represents  the  known  analytical  properties  that  the  series  solution  has  at  large  distances, 
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where  the  series  cannot  simply  be  summed.  In  the  last  section,  the  techniques  used  to  find  the 
full  asymptotic  expansion  are  described. 


2.  Development  of  the  Radial  Equations 


The  equations  for  the  radial  wave  functions  of  the  quadrupole  guide  are  derived  by  making  use 
of  the  angular  symmetry  of  the  quadrupole  magnetic  field.  This  symmetry  is  expressed  as  the 
conservation  of  alignment,  Az  =  Lz  —  Sz  (Golding,  2009;  Hinds,  2001;  Lesanovsky,  2004).  The 
approach  taken  here  is  to  find  eigenfunctions  that  are  common  to  both  the  Hamiltonian,  H  and 
the  alignment,  Az.  This  method  effectively  separates  the  angular  and  radial  dependence  of  the 
problem  leaving  only  a  coupled  system  of  radial  equations  to  solve. 

By  removing  the  angular  dependence,  the  following  pair  of  coupled  radial  equations  for  the 
specific  alignment,  // ;  energy,  s ;  bias  field,  b0 ;  and  field  gradient,  bx ,  can  be  derived: 


8-,K  +  ~d <A-^  +  V2)  =  btPR_, 

P  P 

,  1  (//-1/2)2 
dpR-+~ d  R  -  — — R_+(e-b0)R_=  bxpR+ . 

P  P 


(1) 


The  parameters  and  physical  relevance  of  these  equations  have  been  described  previously  by 
Golding  (2009).  In  the  case  of  a  spin  half  particle  with  integral  orbital  angular  momentum,  the 
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pair  must  be  solved  for  the  allowed  half-integral  values  of  the  alignment,  ju  =  ±  — ,  , 

although  the  solutions  for  negative  alignment  are  simply  derived  using  symmetry  when  the  bias 
field,  b0  is  zero. 

An  important  aspect  of  these  equations  is  that  they  are  coupled.  In  system  1,  R+  and  R  occur  in 

both  equations.  In  other  studies  of  this  system  it  is  assumed  that  the  boundary  conditions  at  the 
origin  can  be  determined  by  using  the  fact  that  as  p  — »  0  the  terms  on  the  right  hand  side  of 

system  1  can  be  ignored  (Hinds,  2000;  Lesanovsky ,2004;  Potvliege,2001).  This  is  equivalent  to 
assuming  that  bx  can  be  set  to  zero  and  limiting  the  results  to  small  radii.  Using  this  assumption 

it  would  appear  that  the  two  equations  reduce  to  an  uncoupled  pair  of  second  order  differential 
equations  for  a  small  region  near  the  origin.  This  is  a  type  of  singular  perturbation  since  the  fonn 
of  the  equations  changes  drastically  when  the  small  parameter  bx  is  set  to  zero. 

Physically,  setting  b{  to  zero  is  like  turning  off  the  trapping  fields.  The  boundary  conditions  for 
the  singular  case  of  bx  =  0  must  be  those  of  a  free  particle  in  which  the  spin  components  are 
completely  uncoupled.  Free  particle  states  are  eigenstates  of  the  spin  operator,  S.  but  they  are 
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not  eigenstates  of  the  full  Hamiltonian.  In  order  to  get  the  proper  boundary  conditions  when  the 
trapping  fields  are  on,  it  is  important  to  first  uncouple  the  equations  and  perform  a  power  series 
analysis  on  the  resulting  fourth  order  differential  equation  using  the  Frobenius  method.  The 
resulting  series  solution  is  evaluated  at  the  origin  and  the  behavior  near  the  regular  singularity  is 
then  handled  correctly  without  making  the  assumption  that  bx  -  0 . 


3.  Uncoupled  Equations 


The  system  of  radial  equations  displayed  in  system  1  contains  an  irregular  singularity  at  infinity. 
In  order  to  develop  an  asymptotic  series  for  large  values  of  the  radius,  the  same  decoupled  fourth 
order  equations  that  are  developed  to  study  the  behavior  near  the  origin  are  used.  By  solving  the 
top  equation  of  system  1  for  R  and  substituting  this  result  into  the  bottom  equation,  the  system  is 

uncoupled  to  produce  the  following  fourth  order  differential  equation  for  R  , 


4  d4 

P  ~r^R+  + 


dp 

2/?4£-f^-  +  2/r 

V  2 


p 


d2 


dp 2 


-/?  + 


-2p\+(5+4f,2  +6fl)p]j-Rt 


(2) 


f  i  \  i  q 

-b2pb  +[s2 -b2)pA  +  b0-2p2s  +  2pb0+-£  p2 -  — p2  -9p-  —  + p 

V  2)2  16 


R=  0 


An  equivalent  equation  for  the  R_  component  can  be  obtained  from  equation  2  by  simply 
changing  the  sign  of  both  p  and  b0  as  well  as  replacing  R+  with  R  .  This  can  be  seen  by 
inspection  of  system  1 . 

Equation  2  is  a  fourth  order  Hamburger  equation  (Ince,  1956)  that  can  be  solved  explicitly  using 
the  Frobenius  series  technique  as  shown  by  Golding  (2009,  2010).  It  has  a  regular  singularity  at 
the  origin  and  an  irregular  singularity  at  oo  .  The  series  approach  is  useful  as  an  expansion  around 
the  regular  singular  point  at  the  origin.  However,  the  Frobenius  technique  cannot  generally  be 
used  to  find  an  expansion  around  the  irregular  singular  point. 

In  equation  2,  the  equation  for  R+  has  been  developed  by  eliminating  R  from  the  coupled 

equations  in  system  1 .  The  behavior  of  the  other  spinor  component  is  obtained  from  the 
following  equation: 
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*(/>)=—■ 

b\P 


d 

d2  „  /  x  J/7 


<3?  yC> 


-^+ (/>)  +  - 


f  1 

//+? 

V  2y 


p 


p 


-  +  (£  +  bo)R+{p) 


(3) 


that  is  easily  obtained  from  the  first  equation  in  system  1 . 

In  order  to  obtain  useful  information  about  the  behavior  of  the  solution  as  p  — »  go  ,  an  asymptotic 
expansion  must  be  developed.  This  requires  first  the  detennination  of  the  controlling  factor  and 
the  leading  behavior  of  the  uncoupled  fourth  order  equations. 


4.  Calculating  the  Controlling  Factor  and  Leading  Behavior 


The  first  step  in  a  solution  of  the  asymptotic  behavior  of  R+  in  equation  2  is  to  calculate  the 

controlling  factor  of  the  expansion.  This  is  done  by  assuming  a  trial  solution  of  the  form, 

R+  =  exp  (.S’  ( pj)  (Bender,  1978).  A  differential  equation  for  .S’(/r)  is  derived  that  must  be 
approximately  solved  for  large  p  .  The  full  equation  is 


f  d  V 

~Ts(p) 

\dp  j 


f  ,  d2  s  „  2/u2  +5/2 V  d 

l  dp2  K  ’  P2  )\dp 


V 


dp' 


f  d2 


\dp 
d4 


s(p) 


p 


-2s  + 


P 


J 


s(p) 


.2  .  d2 


2p  +5/2 
P2 


dp ' 


«W- 


dp 


S  (/?)  +  by  p2  +  bj  -s2  + 


2 sju  -  2 boJu  -1/2 £  —  b0 
P2 


( 


9ju-ju*+  —  - +19/2//2  I p-4  =  0 

v  16  J 


(4) 


By  using  the  several  approximations  that 
simplified  equation  is  obtained  for  S(p) 


d" 


( 


dpn 


S« 


\dpL  j 


as  p  gets  large,  the  following 


(  d  ^ 

j ~S(P 
ydP 


-2s 


d 


-S(p)^  +  b2p2  -  ’  0 


(5) 
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By  treating  equation  5  as  a  quadratic  equation  in  -^-5  (/:?  ),  the  following  four  differential 

dp 

equations  are  derived: 


dS_ 

dp 

dS_ 

dp 


±iy[b^p  +  £. 


(6) 


The  solutions  are 


S(p)  =  ±^-(blP  +  er+C(p) 
s  ^ = ±z  ~£^2 + c  ^ 


(?) 


where  C(p )  is  an  integration  “constant”  to  be  determined.  The  idea  in  the  asymptotic  analysis  is 
that  C(p)  should  be  much  smaller  at  infinity  than  the  first  term  in  S (p)  ■  Using  the  techniques 
described  by  Bender  (1978),  C(p)  is  detennined  by  substituting  the  solutions  in  equation  7  back 


into  equation  4  and  using  the  approximations  that 


d" 


dp" 


C« 


-*-■ C 

dp  j 


to  obtain  the  leading 


behavior.  By  this  technique,  the  factor  C(p)  is  found  to  be  -  —  In  bxp  and  the  leading  behavior  is 
determined  by  the  four  functions 


exp 


2  co 


{bxp  +  s)m 


(vf 


(8) 


in  which  the  factor  co  represent  the  four  roots  of  unity,  ±/',±l .  To  be  consistent  with  the 
asymptotic  nature  of  this  calculation,  the  exponent  is  expanded  and  the  terms  that  become  small 
at  large  p  are  neglected.  The  factor  S  (p)  is  then  given  by 


2co 
3  \ 


(bxp  +  s)m 


2co/3^bl(pf2 


(9) 


This  expansion  of  the  exponent  is  consistent  with  the  normal  solutions  at  irregular  singularities 
discussed  by  Ince  (1956).  However,  this  choice  is  not  unique  and  other  possible  choices  may  or 
may  not  result  in  a  consistent  asymptotic  series. 
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Using  equation  9,  the  leading  behavior  of  the  radial  wave  function  is  given  by  the  four  functions 
that  we  refer  to  as  the  asymptotic  basis  functions 


exp 


(10) 


The  behavior  of  either  spinor  component,  R+  or  R  at  large  p  must  be  asymptotic  to  some 

superposition  of  these  four  asymptotic  basis  functions.  This  is  true  even  though  the  functions  in 
equation  10  were  derived  from  the  equation  for  R+ ,  equation  2.  The  reason  for  this  is  that  the 

equivalent  fourth  order  equation  for  R  only  differs  in  the  signs  of  //  and  b0 ,  neither  of  which 

appear  in  equation  10.  The  leading  factors  for  the  two  components  are  therefore  identical.  The 
actual  dependence  on  //  and  b0  will  reappear  as  more  terms  in  the  asymptotic  expansion  are 

included. 


5.  Stokes  Phenomenon  and  the  Asymptotic  Basis  Functions 


In  the  global  basis,  all  the  physically  allowed  solutions  found  using  the  Frobenius  technique  in 
this  problem  are  of  the  form  (Golding,  2010) 

oo 

/^2Xp!",  (in 

n= 0 


where  a  is  either  an  even  or  odd  integer  that  depends  only  on  the  mode  and  component  of 
interest.  In  the  local  basis,  the  components  are  superpositions  of  similar  series  of  different 
indicial  exponents  <7  and  the  analytic  properties  of  the  combined  series  are  not  as  simple.  For 
this  reason,  the  global  basis  is  used  and  the  series  form  in  equation  1 1  is  taken  as  the  general 
form  for  the  asymptotic  analysis  below. 

Once  this  series  is  extended  to  the  complex  plane  by  letting  p  — »  pe‘° ,  the  coefficients  needed  to 
combine  the  asymptotic  basis  functions  of  equation  10  into  a  proper  approximation  of  the  power 
series  in  equation  1 1  can  be  determined.  This  is  accomplished  by  forcing  both  representations  to 
have  the  same  analytic  properties  for  large  p  .  The  extension  of  the  asymptotic  basis  functions  to 

the  complex  plane  are  also  obtained  by  the  substitution  p  — >  pe"  and  are  given  by 


exp 


L(P’8)  =  - 


2  co 

v3 A 


(Vf 


COS 


V  *  ) 


+  i  sin 


flQ^ 


v  2  JJ) 


(V) 


3/4 


.3(9 


where  a>  =  ±i,  ±  1 . 


(12) 
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The  expression  in  equation  1 1  represents  a  single-valued  function  of  6  in  the  complex  plane.  A 
reasonable  approximation  to  this  complex  function  at  large  values  of  p  can  be  formed  by  using 

various  linear  combinations  of  the  four  asymptotic  basis  functions  in  equation  12  to  create  a 
single-valued  function.  However,  the  behavior  of  the  individual  basis  functions  in  equation  12  is 
dramatically  different  in  different  sectors  of  the  complex  plane.  Different  linear  combinations  of 
the  basis  functions  are  required  in  each  sector  to  approximate  the  series  solutions.  The  different 
representations  in  each  sector  must  in  the  end  be  matched  to  the  adjoining  sectors  so  that  an 
effectively  continuous  function  is  fonned.  This  matching  process  is  handled  using  the  phase 
integral  techniques  designed  to  deal  with  Stokes  Phenomenon  (Bender,  1978;  Heading,  1962; 
White,  2005).  The  phase  integral  calculation  takes  place  on  the  complex  z-plane  shown  in  figure 
1.  The  basic  technique  is  described  in  Heading  (1962)  and  White  (2005)  for  second  order 
differential  equations.  In  figure  1  the  lines  described  as  anti-Stokes  lines  radiate  outwards  along 
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the  lines  6  =  0,  ±  .  These  lines  cut  the  plane  up  into  three  sectors  where  different  forms  of  the 

solution  are  used. 


Figure  1.  Complex  plane  showing  the  positions  of  the  anti-Stokes  lines 
needed  to  perform  the  phase  integral  calculation  that  properly 
combines  the  asymptotic  basis  functions. 

The  anti-Stokes  lines  are  defined  so  that  the  pair  of  exponential  solutions  with  co  =  ±i  in  equation 
12  have  pure  imaginary  phase  and  therefore  do  not  decay  towards  infinity.  Along  another  set  of 
lines  called  the  Stokes  lines  (the  red  lines  in  figure  1)  these  same  oscillatory  solutions  have 
purely  real  exponents  and  one  solution  decays  towards  infinity  while  the  other  grows 
exponentially.  The  real  component  of  the  phase  changes  sign  on  crossing  an  anti-Stokes  line. 
Thus,  an  exponentially  growing  solution  in  one  sector  abruptly  becomes  an  exponentially 
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decaying  solution  in  the  adjacent  sector.  Phase  integral  techniques  are  intended  to  smooth  out 
this  apparent  discontinuity. 

An  important  consideration  in  this  problem  is  that  there  are  four  functions  in  equation  12  that 
must  be  considered.  These  functions  form  two  pairs  and  the  Stokes  lines  of  one  pair  are  the  anti- 
Stokes  lines  of  the  other  pair.  For  example,  along  the  real  axis  where  0  —  0^  the  solutions  with 

co  -  ±i  have  purely  imaginary  phase  and  thus  do  not  decay  towards  infinity  while  the  solutions 
with  co  -  ±1  have  purely  real  phase  and  break  up  into  decaying  and  growing  exponentials. 

The  calculation  proceeds  by  starting  in  sector  7  of  figure  1  with  a  superposition  of  the  form 
+  Bi;\  +  C^\  +  Dq\  ( A,  B,  C,  and  D  are  assumed  to  be  complex  constants).  The 
superscripts  s  and  d  refer  to  the  subdominant  or  dominant  nature  of  the  basis  function  in  the 
sector.  Subdominant  components  decay  towards  infinity  and  dominant  components  grow. 

The  first  step  is  to  determine  a  superposition  in  region  1  that  is  consistent  with  the  assumed  fonn 
in  region  7.  This  is  accomplished  using  arbitrary  Stokes  constants  and  by  adjusting  the 
dominancy  of  the  solutions  as  described  in  both  Heading  (1962)  and  White  (2005).  Since  the 
phase  factor  in  the  denominator  of  the  basis  functions  is  independent  of  p ,  it  can  be  ignored 

except  when  crossing  the  branch  cut.  Once  this  process  is  complete  the  solution  in  region  2  is 
developed  from  that  in  region  using  the  Stokes  constant,  7j .  This  process  proceeds  around  the 

complex  plane  until  after  dealing  with  the  branch  cut,  region  7  is  reached  again.  Referring  to  the 
diagram  in  figure  1  the  detailed  calculation  goes  as  follows: 

1 .  Starting  in  region  7  with  the  assumed  form  A£4i  +  Bg\  +  C^‘t\  +  Dd\ 

2.  Cross  into  region  1  obtaining  A%s+.  +  B£di  +  Cd‘{2  +  (D  +  TtC)d\ 

3 .  Cross  into  region  2  obtaining  ( A  +  T2B ) %s+.  +  Bq1,  +  Cd'\  +  (D  +  TXC ) d‘\ 

4.  Cross  into  region  3  obtaining  (A  +  T2B')2;4i  +  B%s_i  +  (C  +  T2  (O  +  TjC))^  +[D  +  TXC)^X 

5.  Cross  into  region  4  obtaining 

(a  +  T,  B)  C .  +  {B  +  r4  (A  +  T2B))  C  +  (c  +  r3  (D  +  T,C))&  +  (D  +  TtC)  £, 

6.  Cross  into  region  5  obtaining 

(A  +  T1B)&+(B  +  T,(A+T1B))g  + 

(C+T,(D  +  T,C ))  +  ((£>  +  Tfi)  +  Ts  (C  +  T,  (D+ Tfi)))^ 

7.  Cross  into  region  6  obtaining 

((A  +  T,B)  +  Tt  (B  +  Ta(A  +  T,B)))?U  +  (B  +  Ta(A  +  T,B ))£  + 

(C  +  T,(D  +  T1C))Z1  +  ((D  +  T1C)  +  TS(C  +  T3(D  +  T£)))& 

8.  Cross  the  branch  cut  to  get  back  to  region  7  obtaining 

+i({A  +  T2B)  +  T6(B  +  T4{A  +  T2B)))  £,  +  i(B  +  T4(A  +  T2B ))#  + 

i(C  +  T2(D  +  TxC))  +i((D  +  TlC)  +  T5(C  +  T3(D  +  TlC  )))  ^ 
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At  this  point  there  are  two  different  descriptions  of  the  function  in  region  7  that  can  be  made 
equal  by  calculating  consistent  values  of  the  Stokes  constants.  By  this  process  there  is  essentially 
a  different  description  of  the  desired  function  in  each  region  of  the  plane  and  these  functions 
must  all  be  connected  together  so  that  the  patched  representation  of  the  function  in  equation  1 1  is 
as  continuous  and  singled  valued  as  possible.  By  comparing  the  final  expression  obtained  in 
region  7  with  the  expression  that  began  the  calculation  in  region  7,  four  equations  are  obtained 
that  are  sufficient  to  determine  that  all  of  the  Stokes  constants,  T.  are  equal  to  —i . 

Once  the  Stokes  constants  are  evaluated  then  the  behavior  of  the  total  solution  in  the  various 
regions  can  be  used  to  completely  calculate  the  required  relations  between  the  amplitudes, 
A,B,C,D  ,  of  the  asymptotic  basis  functions  needed  to  properly  approximate  the  asymptotic 

behavior  of  the  series  solutions  of  equation  1 1 . 

As  noted  before,  the  series  in  equation  1 1  implies  that  the  extended  function  along  the  negative 
real  axis  must  be  equal  to  eian  times  the  function  along  the  positive  real  axis.  Setting  all  of  the 
Stokes  constants  equal  to  —i  and  ignoring  appropriate  subdominant  terms  that  decay 
exponentially  along  the  real  axis,  the  expression  along  the  negative  real  axis  in  region  1  becomes 
(. A  —  iB )£+i  -  iDc't  +(Z)-iC)4Lj  and  the  similar  expression  along  the  positive  real  axis  is 

.3  n 

d  d  '  •  .  ••  •  ^ 

A%s+i  +  Bg  (  +  C£+1 .  Taking  all  of  this  into  account  and  including  the  phase  shift,  e  4  common  to 

all  of  the  basis  functions  evaluated  along  the  negative  real  axis,  the  following  equation  is 
obtained: 

.(13) 


Equating  the  coefficients  of  common  basis  functions  in  this  expression  the  relations 

.3k 

(A-iB)e"T  =  ei<mC 

-iDe  4  =  eimB  ,  (14) 

.371 

(. D-iC)e’~*  =e‘™A 


A  exp 


-(v)” 

v3i>, 


f 


J 


.3  71 
l~4 


+  Bex p 


2 i  x3/2 


V  3ft, 


+  Cexp 


J 


\ib\  j 


(A-iB)ex  p  — (V) 


3/2 


V3  \ 


-  iD  exp 


J 


-^{bxpf2\{D-iC)exp[^{bxp)m 

\  )  V  J 


are  found.  The  eigenfunctions  described  by  the  series  in  equation  1 1  must  be  finite  and  real 
valued  along  the  positive  real  axis.  This  is  enough  information  to  detennine  that  the  coefficients 
are  given  by 
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(15) 


A  =  ±\A\e'4 

.71 

B  =  ±\A\ e‘~A  . 

C  =  0 

D  =  +  \A\e™ 

The  two  possible  sets  of  solutions  are  necessary  to  account  for  the  sign  of  the  exponentially 
decaying  components  of  the  radial  eigenfunctions  as  the  radial  mode  number  is  increased.  That 
is,  if  the  lowest  energy  eigenfunction  of  a  particular  mode  has  no  zero  crossings  then  D  is  taken 
to  be  positive  allowing  the  bound  component  of  the  wavefunction  to  decay  from  positive  values 
towards  zero.  If  a  is  odd  the  upper  sign  is  taken  and  if  a  is  even  the  lower  sign  must  be  taken. 
The  next  higher  radial  mode  would  require  D  to  be  negative. 


6.  Results  of  Phase  Integral  Calculations 


Using  the  constants  determined  by  the  phase  integral  technique  the  functional  forms  of  the 
leading  behavior  for  each  component  and  each  mode  are  displayed  here: 


R„  =  2- 


cosl  |,/U“  +  Up/ai 


(Vf 


-(-1) 


n+ 5/2  e 


“T  yfcp^+ey/p/bi 


(M 


3/4 


COS 


K  1—2- 


cos 


R  ,  =  -2- 

m2 


“  +£yjp/b i  +  — 

(vf  ‘ 

y  &P*  2  +  e^~pTbx  +  — 


-(-1) 


p+5/2  e 


z4b\PV2+spP+ 


(M" 


3/4 


+  (-l) 


JU+V2  e 


-y!b\Pm+Sy[pibx 


(16) 


COS 


;r 


^2=2- 


j  + \/'P ;  /'i  +  ^ 

(vf 


+  (-!) 


n+n2  e 


-yfbiP3/2+£ylP/b l 


These  four  functions  represent  the  leading  behavior  of  the  spinor  components  and  are  valid  along 
the  positive  real  axis  as  /?  — >  °o  .  The  proportionality  factor  |d  has  been  set  equal  to  one  as  these 

solutions  need  to  be  matched  to  series  solutions  describing  the  behavior  in  the  area  of  the  origin. 
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The  functions,  Rpl  and  Rp2 ,  refer  to  the  spin-up  components  of  type  1  modes  and  type  2  modes 

(Golding,  2009)  and  Rm]  and  Rm2  similarly  refer  to  the  associated  spin-down  components.  The 

phase  choice  in  these  expressions  produces  a  positive  decaying  exponential  for  the  tail  of  the 
bound  component  of  the  ground  state  for  //  =  1  /  2  .  The  overall  phase  is  unimportant  but  the 

relative  phase  between  the  cosine  and  exponential  tenns  in  each  function  is  the  mathematical 
cause  of  the  exponential  decay  normally  seen  in  the  bound  state  components,  R  +  Rm .  A  quick 

check  shows  that  the  sinusoidal  tenns  cancel  out  in  the  bound  components  and  the  exponentials 
cancel  out  when  the  difference  components,  R  -  Rm  are  calculated.  Thus  each  of  the  global 

spinor  components,  Rp ,  Rm  shown  in  equation  12  contains  a  mixture  of  an  unbound  and  a  bound 

component. 

Within  a  given  series  of  radial  eigenfunctions  the  boundary  conditions  are  identical  for  each 
mode  at  the  center  of  the  guide.  As  the  eigenenergies  increase  from  one  radial  mode  to  the  next, 
the  sign  of  the  exponential  factor  must  change  to  account  for  the  increasing  number  of  zero 
crossing  at  higher  energies.  This  can  be  handled  by  multiplying  each  of  the  functions  in  equation 
12  by  a  factor,  -l"1  or  -l"2 ,  the  new  parameters,  nx  and  n2 ,  are  then  the  radial  mode  indices  for 

the  radial  wavefunctions.  These  factors  are  not  included  above  but  can  be  easily  included  when 
working  with  a  full  set  of  modes. 


7.  Asymptotic  Expansion  of  the  Mode  Functions 


In  order  to  obtain  a  more  accurate  form  of  the  mode  functions  than  just  the  leading  behavior 
shown  in  equation  16,  a  full  asymptotic  expansion  must  be  calculated.  This  is  accomplished  by 
solving  equation  2  using  an  assumed  form  that  is  made  up  of  the  product  of  the  leading  behavior 
of  the  mode  times  a  series  consisting  of  inverse  powers  of  p  .  The  coefficients  of  the  series 
expansion  are  then  determined  recursively  by  equating  the  factors  multiplying  each  power  of  p 
in  the  resulting  series  zero. 


Substituting  the  following  expansion  in  equation  2, 

(2  -  - ^ 


K  (p)  =  wi  ( P ) 


cosl-^p 


3/2 


+  £yjp  I  "b\ 


(Vf 


+  w2(p) 


sinl  -  Vv 


3/2 


+  £ 


JpTbi 


(Vf 


(17) 


results  in  a  set  of  coupled  differential  equations  for  w,  and  w2  that  can  be  solved  by  assuming 
that  these  functions  have  the  following  asymptotic  expansions: 

oo  oo 

w,  =  Y.a,p-'2  aniw2  =  Y.c,p-”2 .  (18) 

n= 0  n= 0 
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Solving  the  resulting  coupled  series  expansions  for  the  first  few  coefficients  results  in  the 
following  approximate  fonn  for  w,  when  60  =  0 : 


wl  ~  a0  + 


£2  b 
1/4^— -1/2  0 


2  A 
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yfp 
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and  the  corresponding  form  for  vv2 
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£ 2  l) 

l/4-^r-l/2  ° 
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(20) 


In  this  section,  some  of  longer  expressions  are  written  in  stacked  forms  within  parentheses  to  try 
and  clarify  the  expressions.  Do  not  confuse  these  expressions  with  matrices  as  they  are  to  be  read 
as  a  single  series  of  terms.  In  order  to  maintain  the  leading  behavior  form  found  by  the  phase 
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integral  technique  in  equation  16,  the  coefficients,  a0  and  c0 ,  must  be  chosen  to  produce  the  — 


phase  shift  in  the  cosine  terms  shown  in  equation  16.  This  is  easily  accomplished  by  setting  the 
first  coefficient,  a0  -  -  J=,  and  the  second  coefficient,  c0  =  --J= . 

The  first  few  terms  of  the  asymptotic  expansion  of  R+  are  then 
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cos  |  -y/b -p3l2+eJpTh 
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(21) 


Ignoring  the  second  term  in  each  bracket  as  p  — >  <x> ,  the  expression  reduces  to  the  form  found  in 


n 


equation  16  along  with  the  —  phase  factor  that  comes  from  combining  the  sin  and  the  cos 
tenns.  Similar  expansions  are  found  for  the  other  components  and  the  other  modes  of  the  guide. 
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An  expansion  for  the  purely  exponential  term  shown  in  equation  16  can  also  be  included  by 
straightforward  extension  of  the  trial  function  shown  in  equation  17  to  include  an  exponential 
times  another  asymptotic  series,  vv3  ( pj .  Since  the  exponential  dies  off  so  quickly  this  is  only 

important  for  an  asymptotic  representation  of  the  atomic  wavefunction  relatively  near  to  the 
center  of  the  guide  yet  outside  the  classical  turning  point  of  the  mode  of  interest. 

To  attain  good  accuracy,  a  significant  number  of  terms  are  needed  in  the  asymptotic  expansions 
However,  the  expressions  are  long  and  difficult  to  display  so  only  a  few  terms  are  shown  here. 
The  unbound  component  is  the  component  that  is  pointed  in  the  direction  of  the  local  magnetic 
field.  This  component  sees  the  magnetic  potential  energy  decreasing  towards  infinity  and  is 
repelled  from  the  center  of  the  guide.  The  unbound  component  is  given  by  R+  (/?)  -  (/?) .  It 

decays  fairly  slowly,  is  not  extremely  sensitive  to  the  bias  field,  b0  and  the  first  few  terms  of  its 
expansion  are  displayed  in  equation  22, 

R+-R  ~ 

42  \l244bp -M \42s2 

- 1 - b 

3/4  3/4  i9  4  J  4 

i  P  R  P  (22) 
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- 1 - - - h 
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Fy  is  A 

1  / +  (l /  8 V2 b04  - 1  / 8 42s- X1  + 1  / 3242s4 ) b~ 4  p~114 

bi  V  ' 

The  other  interesting  component  for  any  mode  is  the  bound  component.  This  is  given  by 
i?+(p)  + /?_(/?),  which  decays  exponentially  except  for  a  small  oscillatory  component  in  the  tail 

given  by  equation  23  when  the  bias  field,  b0  is  zero, 

R++R_  ~ 

15  13  17  " 

-l/242jub4,4p^  +l/8V2fWV~T  +  (23) 

15  13  _17 

\/242pbt  ~ll4p  4  + 1  /  S42s2pb1  4  p  4 

Notice  that  this  oscillating  tenn  drops  off  like  — and  is  proportional  to  the  alignment. 

p 
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However,  when  the  bias  field  is  increased  there  is  a  pronounced  effect  on  the  behavior  of  the 
oscillatory  terms  that  make  up  the  bound  component  as  can  be  seen  in  equation  24, 
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(24) 


All  of  the  terms  displayed  in  equation  24  are  proportional  to  the  bias  field,  b0  and  vanish  at  zero 

bias.  Note  also  that  the  radial  dependence  is  much  stronger  in  these  bias  dependent  components. 
These  effects  are  small  at  low  bias  fields  but  are  expected  to  become  important  when  studying 
the  adiabatic  problem  when  it  is  normal  to  include  large  bias  fields  to  eliminate  the  possibility  of 
some  types  of  spin  flips. 


8.  Conclusions 


A  method  for  calculating  complete  asymptotic  expansions  for  the  modes  of  atomic  waveguides 
has  been  demonstrated.  The  technique  can  be  extended  to  high  orders  using  computer  algebra 
techniques.  This  method  complements  the  series  solution  technique  presented  in  Golding  (2010) 
and  produces  an  accurate  representation  of  the  waveguide  solutions  far  away  from  the  guide 
center. 

In  this  system,  some  components  of  the  spinor  wavefunctions  do  not  die  off  exponentially  and 
are  therefore  significant  at  large  distances.  This  is  related  to  the  quasibound  nature  of  the  system. 
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A  good  representation  of  the  nondecaying  components  of  the  eigenstates  at  large  distances  from 
the  guide  center  is  required  so  that  the  calculated  eigenstates  can  be  used  reliably  in  further 
numerical  calculations.  The  asymptotic  expansions  presented  here  provide  this  representation 
and  are  readily  related  to  the  power  series  solutions  developed  in  Golding  (2010).  By  connecting 
those  power  series  solutions  to  the  asymptotic  expansions  developed  here  an  efficient 
representation  of  the  exact  radial  wavefunctions  can  be  obtained.  These  wavefunctions  are 
needed  for  detailed  studies  of  important  properties  of  magnetic  guides  such  as  sensitivity  to  noise 
driven  spin  flips,  importance  of  quantum  Majorana  transitions,  energy  level  dependence  on 
magnetic  field,  and  the  effects  of  guiding  field  imperfections  as  well  as  the  onset  and  departure 
of  adiabatic  behavior. 
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